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■ Abstract 

We propose new construction of the polynomial integrals of motion related 
to the addition theorems. As an example we reconstruct Drach systems and get 
some new two-dimensional superintegrable Stackel systems with third, fifth and 
seventh order integrals of motion. 
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1 Introduction 

The Liouville classical theorem on completely integrable Hamiltonian systems implies 
that almost all points of the manifold M are covered by a system of open toroidal 



C<S ■ domains with the action-angle coordinates / = {ly, . . . , I k ) and u> = (oji, . . . , u n ): 

in 

O; {ij,i k } = {ui,u} k } = o, {i j ,u k } = 5 ij . (l.i; 
oo 

The independent integrals of motion Hi, ... , H n are functions of the independent action 
variables I\ , . . . , I n and the corresponding Jacobian does not equal to zero 

£ : detJ^O, where J i3 - = dHi ^ r " " ' ^ . (1.2) 



dlj 



Let us introduce n functions 



k 

such that 

{tfi,M = £j*( J_1 )jy = V (1.4) 



kj 

fc=l 



If Hamiltonian H = Hi then the (n — 1) functions <p2, ■ ■ ■ , <j>n are integrals of motion 

^- = {H U ^} = 0, j = 2,...,n, 
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which are functionally independent on n functions Hi (/),... , H n (I). 

So, in classical mechanics any completely integrable system is superintegrable system 
in a neighborhood of any regular point of M |14j . It means that the Hamiltonian 
H = Hi has 2(n — 1) integrals of motion H 2 , . . . , H n and fa, . . . , <f> n on any open 
toroidal domain. 

If the action-angle variables are global variables on the whole phase space M and, 
therefore, we have superintegrable systems on M. For instance, the global action-angle 
variables for the open and periodic Toda lattices are discussed in [3]. 

However, in generic case the angle variables u k are multi-valued functions on the 
whole phase space M. If we have k additional single-valued algebraic integrals of 
motion K the trajectories are closed (more generally, they are constrained to an n — k 
dimensional manifold in phase space). 

Any additional integral is a function on the action-angle variables. Since we have 
to understand how to get single-valued additional integrals of motion from the multi- 
valued action-angle variables. In this paper we discuss a possibility to get polynomial 
integrals of motion from the multi-valued angle variables by using simplest addition 
theorem. 



2 The Stackel systems. 

The system associated with the name of Stackel (9j [TO] is a holonomic system on 
the phase space M = R 2n , with the canonical variables q = (qi,...,q n ) and p = 
(pi,---,Pn): 

n 

^ = 7", dpj A dqj , {pj,qk} = 5 jk . (2.1) 

3=1 

The nondegenerate n x n Stackel matrix S, whose j column depends on the coordinate 
qj only, defines n functionally independent integrals of motion 

n 

H k = J2 < < S ~ 1 ^{p 2 3+ U ^))- (2-2) 

3=1 

From this definition one immediately gets the separated relations 

n 

p] = Y,H k S k] -U ] {q J ) (2.3) 

k=l 

and the angle variables 

j= i J Pi j= i J VY2=i H k S kj - Uj{qj) 

It allows reducing solution of the equations of motion to a problem in algebraic geom- 
etry |10| . Namely, let us suppose that there are functions fij and Xj on the canonical 
separated variables 

/'., ".ji'lj)!'.,- ^j = v Mi)i {li,Pj} = hj, (2-4) 
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which allows us to rewrite separated equations (12.31) as equations defining the hyperel- 
liptic curves 

Cj : = Pj(Xj) = Uj(Xj) (y^H.S^) - UjiXM , (2.5) 

where Pj{Xj) are polynomials on Xj. In this case the action variables Ik = H^ (12. 2p 
have the canonical Poisson brackets (11.11) with the angle variables 

«* = £/ -Wf^^^ft). ( 2 - 6 ) 

j=l JAi VPjiXj) j=1 

which are the sums of integrals of the first kind Abelian differentials on the hyper- 
elliptic curves Cj ( 12. 5p [TP} 114] . i.e. they are sums of the multi- valued functions on the 
whole phase space. 

2.1 Addition theorems and algebraically superintegrable sys- 
tems 

In generic case the action variables (I2.6P are the sum of the multi- valued functions dij. 
However, if we are able to apply some addition theorem to the calculation of uj{ (12. 6p 

it 

oJi = $ij{Pj, qj) = Qi(Ki) + const, (2.7) 

3=1 

where G« is a multi- valued function on the algebraic argument K{(p,q) then one will 
get algebraic integrals of motion Ki(p, q) because 

= {H^Q^Ki)} = e' 4 • {Hx,Ki} = 0. 

So, the addition theorems (12.71) could help us to classify algebraically superintegrable 
systems and vice versa. Plane curves with the genus g > 1 are related to elliptic 
and Abelian integrals. Addition theorems of these functions are the content of Abel's 
theorem |15j . 

The main result of this paper is that almost all the known examples of algebraically 
superintegrable systems relate with the one of the simplest addition theorem 

e x e y = e x+y^ Qr ln ( Xi ) + ln ( X2 ) = ln ( XiX ^ ( 2 gj 

associated with the zero-genus hyperelliptic curves Cj 

Cj: /'j I'jiXj) J)\j ■ !lr\, ■ h.r J 1 »• (2.9) 

where fj, gj, hj are linear functions on n integrals of motion Hi, ... , H n . 

In fact, if Sij(X) = 1 then after substitution (12. 9p into (I2.6P one gets the sum of the 
rational functions 

/ / % dX > = V 

J V'ljX, ■ I' 9j 
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or logarithmic functions 



#i = / 1 =dA = /- 1/2 In L,. + 



2 /A + 9 j 



fjXj + gjXj + hj V -v'7/ 

For the hyperelliptic curves of higher genus one gets elliptic functions, which have more 
complicated addition law [15]. 

In order to use addition low (12.81) we have to make the following steps: 

• We have to apply canonical transformation of the time that reduce n-th row of 
the Stackel matrix to the canonical Brill-Noether form [TTJ [T5] 



Snj — 1 j 3 — 1 > ■ • • n 

such that 

CO 



•n = — V/ = dA, (2.10) 



N 

3=1 

where N is normalization, which is restored from {u> n , H n } = 1. 

• Then we have to use addition law (I2.8P for construction of the polynomial in 
momenta integrals of motion 

• We have to make inverse transformation of the time, which preserves the polyno- 
mial form of integrals such us it depends on q variables only [HI [13] . 

Let us consider construction of polynomial in the momenta integrals of motion at 
n = 2. 

2.1.1 Case f x = f 2 = 

If fi,2 = we have 

1 a pto) 1 Pl Ul P2 U 2 K 

uj 2 = -y I — : dA = 1 = , (2.11) 

A J^[J V9j X + h j 91 92 9192 

where 

K = g 2 piu x + gip 2 u 2 

is the polynomial in the momenta integral of motion of the first or the third order. 
Remind that g\^ 2 are linear functions on the Stackel integrals i?i,2, which are the second 
order polynomials on px,2- 

Example 1 Let us consider the two-dimensional Stackel system defined by two Rie- 
mann surfaces 

Ci, 2 : /i 2 = ^1,2 (A) = (H, ± H 2 )X + a h2 , 
and substitutions (12.41) 

/'.; 'IjI'j- x i '/; • 
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which give rise to the following separated equations 

p\ 2 = H 1 ±H 2 + ^. 

<?1,2 

Integrals of motion H\ and H 2 are solutions of these separated equations 

pi ± pi a x a 2 



Hi: 



2 2qf T 2qV 



(2.12) 



which coincide with integrals of motion for the two-particle Calogero system after an 
obvious point transformation 

X= , Px=Pl~P2, V = = , Py=Pl+P2- (2.13) 



The angle variables (12. lip read as 



i r* 



^1,2 



l 

,dA±- 



47 ^[E x + H 2 )X + «i 47 y/fHl - H 2 )X + a 2 

- H 2 ) TP2q2(H 1 + H 2 ) 
2(H 1 + H 2 )(H 1 -H 2 ) 

The corresponding cubic integral of motion K (12.111) is equal to 

K = 2{H l + H 2 ){H X - H 2 )u 2 = -pMHi - H 2 )+p 2 q 2 (H x + H 2 ). 
The bracket 

{K,H 2 } = 2(Hl-H!) 
is easily restored from the canonical brackets (II. ip . 

2.1.2 Case f x = k\f and f 2 = k\f with integer k Xj2 
In this case, we have 



dA 



0J 2 



1 2 

3=1 
1 



2fc 1 A; 2 v7 



In 



y/fj\ 2 + gj\ + h 



2fc,-v7 



where 



So function 



2fc lv ^/ 



p x u x 



P' 

r 2 



+ p 2 u 2 



dX 



2k 2 ^f 
2 fj v i(0j) + 9j ■ 



\=Vj{qj) 



is the generating function of polynomial integrals of motion K m and K% of the m-th 
and t-th. order, respectively. The values of m and t depend on values of kip- 
As above, the algebra of integrals ifi j2 and K m may be restored from the canonical 
brackets (Ojl . 

Example 2 Let us consider the two-dimensional Stackel system defined by the two 
Riemann surfaces 

Cx,2 : /i?, 2 = P1AK2) = kl 2 \\ 2 + /5 li2 A li2 + (Hi ± H 2 ) 

and substitutions (12.41) 

H = Pj, x i = Qj , 
which give rise to the following separated equations 

Pl,2 = k\ 2 q\ 2 + (3i, 2 qi, 2 + (Hi ± H 2 ) . 

Integrals of motion Hi and H 2 are solutions of these separated equations 

_ p\ ± v\ kfqj ± k\q\ (5 x q x ± (3 2 q 2 
Hi,2 - — 2 2 ' (2 ' 15) 

which coincide with integrals of motion for the harmonic oscillator. The same Stackel 
system coincides with the Kepler problem after well-known canonical transformation of 
the time, which changes the row of the Stackel matrix [IT] . 
The angle variable (12.1 ip reads as 

u 2 = -- [ r , 1 =dA- r , 1 =dA 

2 \J ^kj\ 2 + (3i\ + (Hi + H 2 ) J ^/k 2 2 X 2 + (3 2 X + (Hi - H 2 ) 



2 yfcf 



Using suitable branches of J k\ one gets function (12. 14j) 



, + ( t 2k 2 2 q 2 + P 2 \ kl 

$ = exp(2A; 1 A;2u; 2 ) = [px p 2 ' 



2ki J V. 2A; 2 

which generates polynomial integrals of motion K m and Kg. As an example, if hi — 1 
and k 2 = 3 one gets third and fourth order integrals of motion K m and K% respectively. 

Of course, oscillator is one of the well-studied superintegrable systems, which was 
earmarked to illustrate generic construction only. Additional theorem (I2.8P allows us 
to get a huge family of the n-dimensional superintegrable systems, which have to be 
classified and studied. 
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2.2 Classification 



In order to classify superintegrable systems associated with the addition theorem (12. 8p 
we have to start with a pair of the Riemann surfaces 



C r . t i 2 = P j {\) = f\ 2 + g j \ + h j , J = 1,2, 



(2.16) 



where 



f = aH l + f3H 2 + 1 , g J =a 9 ] H l + l3 g J H 2 + 1 g j , hj = a)E x + %H 2 + 7 j\ 

and a, j3 and 7 are real or complex numbers. 

In order to use addition theorem (12. 8p we fix last row of the Stiickel matrix. Namely, 
substituting 

S 2j (X) = Kj (2.17) 

into the (I2.3H2.6I) one gets 



r) 



2.1 



S 2j { qj )d qj 



P., 



,/p~ 



K 3 f- 1 ' 2 In U + 



2x/7 



(2.18) 



so the angle variable 
1 



ijj 2 = —= In 



PlUl + 



Pj \ K1 



p 2 u 2 + 



PJ 



dp 3 {\) 



dX 



\=Vj(qj) 



is the multi-valued function on the desired algebraic argument 



p 2 u 2 + 



r 2 



2V7 



If Ki j2 are positive integer, then 



K 



2Vf 



K1+K2 



(2.19) 



is the generating function of polynomial integrals of motion K m and Kg of the m-th 
and m ± 1-th order in the momenta, respectively. 
As an example we have 



K m = 2{ Pl u 1 P 2 '+p 2 u 2 P[) 
K m = 2P' 2 (2p 2 u 2 P[ + p lUl P^j + 8p lUl p 2 u 2 f 



Ki = 1, K 2 = 1, 
Ki = 1, K 2 = 2, 

K m = 2P' 2 (Zp 2 u 2 P[ + pmxP^ + 8p 2 u 2 (p 2 u 2 P[ + 3 PlUl pf) f, Kl = 1, k 2 = 3, 

(2.20) 

where m = 1,3, m = 3,5 m = 3, 7, because P[ 2 and / are linear functions on H\ 2 , 
which are the second order polynomials on momenta. 
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The corresponding expressions for the K% look like 



K e = P[P^ + A Pl p 2Ul u 2 f, 

K t = P[P' 2 + Afp 2 u 2 (p 2 u 2 P[ + 2 Pl u l P! 2 ) , 

y d/3 



Kl = 1, « 2 = 1, 
Kl = 1, K 2 = 2, 

^ = P^Pf + 12P^ 2 « 2 (P 2 'p lWl + P^u,)/ + 16 Pl u lP 3 2 u 3 2 f 2 . K\ = 1, k 2 = 3. 



(2.21) 

The imposed condition (12.171) leads to some restrictions on the functions Vj(qj) and 
Uj(qj). In fact, substituting canonical variables (12.41) into the equations (12 . 1 6[) we obtain 
the following expression for the Stackel matrix 



/ av I + afvi + a\ av 2 + a 2 v 2 + a 2 \ 



S 



U 7 



L 1 U 2 

(3v 2 + pf Vl + ft (3v 2 2 + (3 9 2 v 2 + (3% 



V 



det S^O. 



(2.22) 



So, for a given « 1>2 expressions for # 2 j (12.181) yield two differential equations on functions 
u, v and parameters (3: 



S 2j { qj 



Kjiijc'j = pvj + rjrj + S'j. j = 1,2. 



(2.23) 



I 


(3 = 0, 


^ = 0, 




= Qj, 




1 


II 


/? = o, 


# = 0, 




= 1, 






III 


(3 = 0, 


/? = 0, 


w j 


= 1, 







For the Stackel systems with rational or trigonometric metrics, we have to solve these 
equations in the space of the truncated Laurent or Fourier polynomials, respectively. 

Proposition 1 If Kj ^ equations $2.23\) have the following three monomial solutions 



(2.24) 



up to canonical transformations. The fourth solution ( IV) is the combination of the 
first and third solutions for the different j 's. 

In order to prove this fact we can substitute u = aq m and v = bq k into the (12.231) and 
divide resulting equation on q m + k - 1 

abkn = b 2 (3q k+1 ~ m + q l - m b$ 9 + q 1 '^^. 

Finally, we differentiate it by q and multiply on q m 

= -( m -k- l)b 2 f3q k - (m - l)b(3 9 -{m + k- l)(3 h q~ k . 

Such as b ^ and k ^ one gets three solutions (12.241) only. 
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Then we suppose that after some point transformation 

x = z x (q), y = z 2 (q), , 2 

Px = wu{q)Pi + wi 2 (g)j?2, P y = w 2 i(g)pi + w 22 (q)p 2 , ' 

where w^- 7^ 0, kinetic part of the Hamilton function Hi = T + V has a special form 

T = [S^ijP] = Su(x, y)pl + gi2(x, y)p x p y + g22(x,y)pl , 
where g is a metric on a configurational manifold. For instance, if we suppose that 

T = J2(S- 1 ) ljP * = PxPy , 
then one gets the following algebraic equations 

W11W21 = [S~ ) n , W12W21 + W11W22 = 0, W12W22 = (5' _1 ) 12 (2.26) 
and the partial differential equations 

{x,Px} = {y,p y } = 1, {Px,y} = {p y ,x} = {p x ,p y } = 0. (2.27) 

on parameters a and functions Zi )2 (qi,q 2 ), Wkj{qi, 92)- 

The remaining free parameters determine the corresponding potential part 

of the Hamiltonian V(x, y). In fact, since integrals Hi i2 is defined up to the trivial shifts 
Hk — > Hk + Cfc, our potential V(x,y) depends on three arbitrary parameters only. 

Summing up, in order to get all the superintegrable systems on a complex Eu- 
clidean space E 2 (C) associated with the addition theorem (12.81) we have to solve equa- 
tions ( I2.23f2.26f2.27l) with respect to functions Uj(qj), Vj(qj), z Xy2 {q x , q 2 ), wjy(5i, 92) and 
parameters a and (3. 

Example 3 Let us consider second solution from the list (12.241) at k x = 1 and k 2 = 2. 
In this case equations (12. 26112. 27p have the following partial solution 

, 1 0\ x = y/qiq 2 , p x = ^Pi + ^p 2 , 

7, 7 2 J y = p y = ^q\q2Px - 2 ±=P2 • 



q 2 ' vhihzi-l 2 ^, 

Adding potential terms one gets two Riemann surfaces 

Ci : [i 2 = H 2 \ 2 - H x \ - 71, and C 2 : fi 2 = H 2 \ 2 - 2 72 A + 4 73 

where fij = Pj, \j = qj 1 . Solutions of the corresponding separated variables are the 
second order Stakel integrals of motion, which in physical variables look like 

rr , , 72 , 73 rj (PxX-p y y) 2 [X X 

Hi =p x p y + Kxy + —== + — , H 2 = — y - 72W--73-- (2.28) 

vxy 3 y 2 4 v v y 

It is new integrable systems, which is missed in the known lists of superintegrable 
systems [DEI HE]. 
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In this case integrals of motion K 5 (I2.20p and Kg (12.211) are the fifth and sixth order 
polynomials in the momenta, respectively. Of course, we can try to get quartic, cubic 
and quadratic integrals of motion K 4 , K 3 and K 2 from the recurrence relations 

K 5 = {K 4} H 2 }, K 4 = {K 3 , H 2 }, K 3 = {K 2 , H 2 } (2.29) 

and the equation {Hi, Kj} = 0, j = 4, 3, 2. 

Solving first recurrence equation one gets quartic integral of motion 

K 4 = -A(xp x - yp y ) (p x + jiy ) + —-^ — + 



xy y/xy y 

The other recurrence relations (I2.29P have not polynomial solutions. Of course, we can 
try modify recurrence relations 

{K 3 , H 2 } = K 4 + F4OH1, H 2 ), {K 2 , H 2 } = K 3 + F 3 (Hi, H 2 ) 

in order to get cubic and quadratic integrals of motion. However, using ansatz 

K 2 = h 1 (x 1 y)p 2 x + h 2 (x,y)p x p y + h 3 (x,y)p 2 y + h 4 (x,y) 

we can directly prove that there is not the additional quadratic integral of motion, 
which commute with Hi (I2.28[) . 



3 The Drach systems 

In 1935 Jules Drach published two articles on the Hamiltonian systems with the third 
order integrals of motion on a complex Euclidean space E 2 (C) with the following Hamil- 
ton function [2] 

H 1 =p x p y + U(x,y). (3.1) 
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Up to canonical transformations x — > ax and y — > by the corresponding potentials look 
like [3 [12]: 

(a) U = — + f3x ri y r2 + ^x r2 y ri , where r 2 + 3r 7 - + 3 = , 

xy j ■> 

(b) u = -?=+ , g + 7(2/ + x) 



xy (y — x) 2 y/xy (y — x) 2 



TT a g 7 



(y — x) 2 (y + x) 2 ' 
(d) ^ = , ° + : g + 7 " 



1) + 1) Vx 2 - l 2 



^ o 2x 2 + 1 7X 

(/) 17 = ax2/ + /%^== + 



v/x^+T Vx 2 + 1 
(#) U = 7 — ; — v? + PKV ~ x ) + 



(h) U = (y+^)-^ 3 



2 lAmxy i m 2 x 2 

- I ( i 1 ! \ 1 1 — ///.//.»! i // 



+ x) 2 3 

a + /3 (y — mx/3) + 7 (y 2 + 

(Jfe) C/ = mT 1/2 + /3xjr 1/2 + 7 x , 
(0 ^ = a( 2/ -^)+/3x- 1 / 2 + 7x- 1 / 2 (i/-px). 

Such as Drach made some assumptions on the form of the third order integrals of 
motion K$ ( 13.61) in the calculation it is not immediately clear whether the obtained list 
is complete. 



3.1 Non-separable systems. 

Let us discuss the Drach systems, which can not be reduced to the Stackel systems by 
any point transformation of variables. 

The first system (a) is non-Stackel system related to the three-particle periodic Toda 
lattice in the center-of-mass frame [12J and there are global action-angle variables [3]. 

The (h) system is reduced to the Stackel system by non-point canonical transforma- 
tion and, therefore, existence of the third order integral of motion is related with this 
non-point transformation [12]. Later this system has been rediscovered by Holt ^j. 

For the (k) case in the Drach papers [2] we can find Hamiltonian 

H[ k) = p xPy + ay~ 3/2 + f3xy' 3/2 + 7X (3.2) 
and the following cubic integral of motion 

rAk) 1 \ ( ^77 dH\ 

Ky = 6w(x,y) l—p y -p x —\-P(p x ,py,x,y), (3.3) 



11 



where 

p = 3p 2 x p y , w = -y. 

It is easy to prove that {H[ k \ ^ and, therefore, we have to suggest the possibility 

of a small mistake in the Drach papers [2]. 

Following to [12], if we solve equation {p x p y + U(x,y), K^} = with respect to 
U(x, y), then one gets our case (k) 

Hi = PxPy + ay~ l/2 + (3xy~ 1/2 + jx. 

On the other hand, we have proven directly that the Hamiltonian H[ (13.21) has one 
second order integral of motion 

=p 2 x - + 2iy (3 . 4) 

and has not cubic integral of motion. Moreover, it is easy to see that quadratic inte- 
grals of motion #f 2 } fj3.2U3.4j) can not be reduced to the Stackel integrals by any point 
transformation of variables. 

Below we do not consider (a) and (h) systems and consider (k) case in our notation 
only. 



3.2 Classification. 

For the Drach systems Ki = k 2 = ±1,1/2 and 

$ 2 = exp(2 v / 7^ 2 ) = {l^JJ +PlU ^)(jJJ + P 2U2 

= L( Ki +y/jK m ) (3.5) 

may be considered as the generating function of the polynomial integrals of motion 
(I2.20ti2.2in 

K m = 2 (pxui P' 2 + p 2 u 2 P{j , m=l,3, (3.6) 

K e = PiP^ + A PlP2 u lU2 f, £ = 2,4, 

of the m-th and £-th order, respectively. It's clear that m = 1, 3 and i = 2, 4, because 
P[ 2 and / are linear functions on H\ 2 , which are second order polynomials on momenta. 

We have to underline, that we use different Ki = k 2 2 = ±1, 1/2 for the agreement 
of K m (13.61) with the initial Drach integrals of motion [2] only. 

One gets third order polynomial integral of motion K 3 (13. 6p if and only if P[(X) or 
P 2 {\) depends on Hi or H 2 . It leads to the additional restrictions on a's and /3's 

tS /0 for ^ 01 3 ^ (37) 

fe=i K 

In order to get all the superintegrable systems on a complex Euclidean space ^(C) 
associated with addition theorem (12. 8p we have to solve equations (I2.23|2.26|2.27l) and 
(13. 7p at Ki j2 = 1 with respect to the functions Uj(qj), Vj(qj), Zi j2 (qi,q 2 ), v?kj(qi,q2) and 
parameters a and (5. 
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Proposition 2 The Drach list of the Stdckel systems with the cubic integral of motion 
Ii3. 6]) associated with the addition theorem A2.8\) is complete up to canonical transfor- 
mations of the extended phase space. 



The results of corresponding calculations may be joined into the table: 





Ci o 112.160 


subs. UM 


21 9 (I2.25P 


s 


b 


/i 2 = H 1 X 2 + {H 2 + 2a)X - (3 + 2 7 
fi 2 = H X X 2 + (H 2 - 2a)X -0-27 


fij = pjqj 


„ _ (gi± 92 ) 2 
-^1,2 — 4 


(f'l) 


c 


„^A> + (V§)a-/J 


^ = 
S' = 1 


_ 91 ±92 
Z h2 - 2 


(f7) 


d 


^2 = # 2 A 2 - v/8(a + /3)A + E x - 2 7 
/x 2 = # 2 A 2 - ^(a - /3)A + E x + 2 7 


/ij = p,j 


~ _ ll+ll 
Zl 2Q1Q2 

z 2 = qiq 2 




f 


,^^_(2_ ffl) A-^f 

a 2 - Q + ^i) A-^ + ^ 
^ V2 V 4 2 


Xj = qj 1 
Kj = -1 


. _ 91-92 
Zl 2 v ^Tq 2 - 

^2 = y/qiQ2 


f f ? ^ 


e 


/i 2 = HxX 2 + 2(/3 + 7 )A + H 2 + 2a 
fi 2 = HxX 2 - 2(j3 - 7 )A + H 2 -2a 


fij = Pj 


„ _ (9l±92) 2 
^1,2 — 4 


(??) 


k 


ft 2 = JA 2 + (P + A + H 2 + a 
2 

/i 2 = ^A 2 + (/? - ifi)A + H 2 -a 

— 


Aj = qj 

Kj = 1 


.* _ 91-92 

Zi - 2 

_ _ (91+92) 2 
^2 - 4 


(91-92) 


g 


/i 2 = -|A 2 + (^ + #2) A + a 
p 3 4 4 8 


A*i=Pi<3i, Al=9i 
H2=2pi, A 2 =(?2 

re 3 = 2 


_ 91 ±92 
Z h2 - 2 


(f7) 



Similar to the oscillator and the Kepler problem, the Kepler change of the time t — > t, 
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where 

dt = v(q)dt, v(q) = —^, (3.8) 

det S 

relates the Drach systems (b),(d) and (e) with the systems (c),(f) and (k), respectively. 
Here S are the Stackel matrices for (b),(d) and (e) systems and S are the Stackel 
matrices for (c),(f) and (k) systems. Remind, that these matrices S and S have different 
first row only, see [TTJ | 



3.3 Integrals of motion 

Using definitions (I3.6P we can prove that integral of motion K± (I3.6P is the function on 
Hi , H 2 and K 3 

K\ = lGhhf 2 + (Kl - Ah l9 2 2 - Ah 2 gl)f + gfgl (3.9) 

Substituting this expression into the definition of $2 (13.51) we can get integral K 3 as 
function on the action-angle variables Ii j2 = #1,2 and u 2 . 

As usual, polynomial algebra of integrals of motion and K 3 follows from the 
canonical brackets (II. ip 

{H U H 2 } = {H 1 ,K 3 } = {H 1 ,K 4 } = 0, 

{H 2 ,K 3 } = 8K A , {H 2 ,K 4 } = 5fK 3 , (3.10) 
{K 3 ,K 4 } = F Z (H U H 2 ,K 3 ), 

where 

• 5 = 4, Fj = lQf{gih 2 + g 2 h x ) - ^g x g 2 {gx + g 2 ) for b,c cases; 

• 5 = -2, F H = K 3 + 32fhih 2 - A(g\h 2 + g\h x ) for d,f cases; 

• 5= 2, Fjjj = 4f(gf + gj) - + h 2 ) for e,k cases; 

• 5= 4, F IV = 2f(gf + 8gih 2 ) -8f 2 h - Ag x gl for g case. 

The difference in the values of 5 is related with the difference in k's which has been 
defined by the Drach integrals of motion [2]. Of course, we can put K\ = k 2 = 1 in all 
the cases such that 5 = 2 and polynomials F look as in (14.11) . It reduces polynomials 
P\ t 2 in the table only. 

As above (I2.29p . we can try to find another second order polynomial integral of 
motion K 2 from the equations 

K 3 = {H 2 ,K 2 }, {H 1 ,K 2 } = 0. 

Solutions of these equations 

K 2 = ±(2p 1 p 2 u 1 u 2 + 2^if 2 / + v x g 2 + v 2 g^j = ± Ki 2 j? 192 (3.11) 
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have been found in [7] in framework of the Lagrangian formalism. The algebras of 
quadratic integrals of motion Hi 2 and K2 have been considered in pQ. 

There is some opinion that all superintegrable systems with quadratic (linear) inte- 
grals of motion are multiseparable, i.e allows the separation of variables in the Hamilton- 
Jacobi equation in at least two different coordinate systems on the configuration space 
[B|. The some of the Drach systems may be considered as counterexamples associated 
with the Lie surfaces [I]. 

Namely, three of the superintegrable Drach systems with quadratic integrals Hi,H%, K2 
are separable in the one coordinate system on the configuration space only. 

Proposition 3 For the (b) and (c) Drach systems integrals of motion Hi,K2 are sep- 
arable in the coordinates 

x = — , V = qiq 2 , 
2?i 

and the corresponding separated relations do not allows us to get cubic integrals of 
motion. 

For the (e) and (g) cases integrals of motion Hi, K 2 are separable in the coordinates 

Qi ~ <?2 (qi + <?2) 2 

* = -2-. y = — 4 

and we can use the corresponding separated relations to the construction of the cubic 
integrals of motion. 

For the (d), (f) and (k) cases quadratic integrals of motion Hi and K2 do not 
separable by the point transformations. 

We can prove this Proposition by using computer program from [5] . 



3.4 The (1) system. 

Without lost of generality we can put p = — 3 in the (1) case. Substituting this Hamil- 
tonian into the computer program from [5] one gets the separated variables 

(9i - Q2) 2 (qi + <? 2 ) 2 



x — — ^ — > V 

and the corresponding separated relations 

p) = Pjfe) = -4aq* =F 8^2^ + m x q\ T 4^/% + H 2 , J = I, 2. (3.12) 

which give rise to one hyperelliptic curve = P(X) at fi — Pj and A = ig^. 
The angle variable 

1 T 1 dA 1 H 2 dA 

uj 2 = 



2J V /P(A) 2J ^PJZxj 

i r qi dA 1 r q2 dA 

2j VTO~2i 



(3.13) 
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is a sum of the incomplete elliptic integrals of the first kind on the common hyperelliptic 
curve. According to [T5l [2] there is addition theorem and additional cubic integral of 
motion 

K 3 = 2(P{p 2 + P 2 'p 1 ), 
which looks like as the Drach integral (13.61) . but in this case functions 

, 2 d P{±q 1>2 ) 



PL = (<?1 + Q2) 



dqi,2 (qi + <?2) 4 



have completely another algebro-geometric explanation. All the details will be pub- 
lished in the forthcoming publications. 

As sequence, the algebra of integrals of motion H\p and K 3 differs from the corre- 
sponding algebras for other Drach systems related with another addition theorem. As 
an example, the recurrence chain Kj + \ = {H 2 , Kj} terminates on the fourths step only 

K 7 = {H 2 , K 6 } = -A80K 4 K 3 + 256H 2 K 3 - 768aH 2 K 3 - 3072(3-fK 3 . 

The solution of the inverse recurrence chain looks like 

K 2 = v\ + 2ax - A^yfx. 

It is interesting, that the algebra of quadratic integrals of motion H\ i and K 2 is one of 
the standard cubic algebras jpQ and we do not explanation of this fact. 

4 New superintegrable systems on zero-genus hy- 
perelliptic curves at K\ = 1 and /C2 = 2, 3 



Let us put Kx = 1 and k 2 = 2 in ( 12. 23ft and try to solve equations (12.26p - fl2.2Tp . Here 
is one superintegrable system with cubic additional integral K m ( 12.201) and quadratic 
integral K t (12T2T|) 



Vni = 71 (3z + y)(x + 3y) + >y 2 (x + y) + 73(2 - y) 
and seven systems with the real potentials 



V^ = 7 1 (3x + y)(a; + 3|/) + 



72 



V, 



(i) 



ii 



V, 



(2) 



lixy + 
7i 



72 



'x 3 y 



+ 



73 

X 2 



(x + yf 



+ 



7a 



xy ar x J 



(3) 



7i 



n 



+ 



72 



xy y^y 



+ 



73 



Vi 
V, 



(1) 

11 

(1) 

i\ 

(2) 
IV 



, 722/ . 732/ 
lixy + —3- + —r 

X X 



2-5/4^3/4 
3 



72 



y/y) 



y/xy ^Jxy 
-yxxy + 72(2 - y) + 



+ 

7s ■ 



73 



xy(y/x + y/y)' 



x + yf 



S 



S 
S 

s 
s 

s 
s 
s 



a b 

.11, 



a/qi b/q 2 
1/31 2/92 

1/31 1 
1/9? V32 2 

1 
1/3? 4/gf 

1 

1/3? 4/g| 

4/g 2 \ 
1/9? ihl ) 

1 q2/4 
1/91 1 



a/gi b 
I/91 1 



16 



for which integrals of motion K m and Ki (I2.201|2T2~1~1) are fifth and sixth order polyno- 
mials in the momenta. 

Solution of the equations K m = ±{H 2 , K m _i} and {Hi, K m ^i} = looks like 

#m-i = / LpiP2U 1 u 2 (2fv 2 + g 2 ) + Av 2 {2fvi + gi) (fv 2 + g 2 ) + {Afh 2 + g\)v x 
= 4/i 2 (^ + P[) - (Afh 2 - g 2 2 )Xi - Ah 29l . 

It is additional integral of motion, which is second order polynomial for the system with 
potential Vm and fourth order polynomial in the momenta for the other systems. 

Now we present some superintegrable Stackel systems at Ki — 1 and k 2 = 3. Here 
is one system with cubic additional integral K m (12.201) 



Vm = 7i (z + 2y)(2x + y) + l2 (x + 2y) + l3 (2x + y), S = ( f \ 
and seven systems with the real potentials 



iiy> 



_ . 72 73 c _ ( 1/gi 



Vr/ - 71^2/ + x 2/3 y 4/3 + x l/3 y 5/3 ' 5 ~ V 9/«2 . ' ' 

,(2) _ 7i . l2s/y 732/" 



2 



o 1 



v {3) _ _7l_ 72 73 Q = ( 1 

V ^ i " x 4/3 y 2/3 " l " x 7/6^5/6' ° I V«i 9 /</2 j ' 



Vg) = 7i(^ 2 - 5x^ + 4^/) + ^ + ^ , 5= ( 2 f 2 f ) , 



73 q / a/gi b 



V IV = lx {x + hy)(hx + y)+ l2 (x - y) + , S = \ 

for which integrals of motion K m and Kp fl2.20H2.2Tl) are seventh and eights order 
polynomials in the momenta. 

The algebra of integrals of motion H\ 2 and K m (I2.20p is the fifth or seventh order 
polynomial algebra in terms of the coefficients of the hyperelliptic curves 

{H 2 , K m } = 2K e , {H 2 , Ke} = 2fK m , {K m , K t } = ±F Z , 

where polynomial Fz depends on the type of solution (12.241) only: 





= 2(Afh 2 


-gi 


F H 


= A(Afh 2 


-gi 


Fin 


= 4(4A 2 


-gi 


Fiv 


= 2(4//i2 


-gi 



^2 



(4.1] 



Here choice of sign + or — depends on k's. 
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As above, the Stackel transformations (I3.8P relate systems associated with one type 
of the solutions (12.241) . whereas algebra of integrals of motion is invariant with respect 
to such transformations. 

The complete classification of such superintegrable systems requires further investi- 
gations. 



5 Conclusion 

We discuss an application of the addition theorem to construction of algebraic integrals 
of motion from the multi-valued action-angle variables. 

We propose new algorithm to construction of the superintegrable Stackel systems 
associated with zero-genus hyperelliptic curves. It allows us to prove that the Drach 
classification of the Stackel systems with cubic integral of motion (13.61) associated with 
the addition theorem (I2.8P is complete. Moreover, we present some new two-dimensional 
superintegrable systems with third, fifth and seventh order integrals of motion. 

The proposed method may be applied to construction of the higher order additional 
polynomial integrals of motion for the n-dimensional superintegrable Stackel systems 
on the different manifolds. On the other hand, we prove that there are some superin- 
tegrable systems, which miss out of this construction. It will be interesting to study a 
mathematical mechanism of the appearance such superintegrable systems. 

The research was partially supported by the RFBR grant 06-01-00140. 
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